Non-critical squeezing in 2-transverse-mode optical parametric oscillators 
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Abstract 

In this article we explore the quantum properties of a degenerate optical parametric oscillator when it is tuned to the first 
family of transverse modes at the down converted frequency. Recently we found [Phys. Rev. Lett. 100, 203601 (2008)] that 
above threshold a TEMio mode following a random rotation in the transverse plane emerges in this system (we denote it as bright 
mode), breaking thus its rotational invariance. Then, owed to the undetermination of the mode orientation, we showed that the 
phase quadrature of the transverse mode orthogonal to this one (denoted as dark mode) is perfectly squeezed at any pump level 
and without an increasing of the fluctuations on its amplitude quadrature (which seems to contradict the uncertainty principle). 
In this article we go further in the study of this system and analyze some important features not considered previously. First 
we show that the apparent violation of the uncertainty principle is just that, apparent, as the conjugate pair of the squeezed 
quadrature is not another quadrature but the orientation of the bright mode (which is completely undetermined in the long 
term). We also study an homodyne scheme in which the local oscillator is not perfectly matched to the dark mode, as this 
could be impossible in real experiments due to the random rotation of the mode, showing that even in this case large levels 
of noise reduction can be obtained (also including the experimentally unavoidable phase fluctuations). Finally we show that 
neither the adiabatic elimination of the pump variables nor the linearization of the quantum equations are responsible for the 
remarkable properties of the dark mode (things which we prove analytically and through numerical simulations respectively), 
which were simplifying assumptions used in [Phys. Rev. Lett. 100, 203601 (2008)]. These studies show that the production of 
non-critically squeezed light through the spontaneous rotational symmetry breaking is a robust phenomenon. 



PACS numbers: 42.50.Dv, 42.50.Lc, 42.50.Tx, 42.65.Yj 

I. INTRODUCTION 

One of the most amazing predictions offered by the 
quantum theory of light is what has been called vacuum 
fluctuations: Even in absence of photons (vacuum) the 
value of the fluctuations of some observables are different 
from zero. These fluctuations cannot be removed by im- 
proving the experimental instrumental, and hence, they 
are a source of non-technical noise {quantum noise) which 
seems to establish a limit for the precision of experiments 
involving light. 

During the late seventies and mid-eighties of the past 
century ways for overcoming this fundamental limit were 
predicted and experimentally demonstrated pj. In the 
case of the quadratures of light (equivalent to the position 
and momentum of a harmonic oscillator), the trick was 
to eliminate {squeeze) quantum noise from one quadra- 
ture at the expense of increasing the noise of its canon- 
ically conjugated one, in order to preserve their Heisen- 
berg uncertainty relation. States with this property are 
called squeezed states, and they can be generated by 
means of nonlinear optical processes. In particular, the 
largest levels of squeezing are obtained by using nonlin- 
ear resonators operating near threshold, this level being 
degraded quickly as one moves away from that critical 
point. 

Up to date, the best squeezing ever achieved is a 
90% of noise reduction with respect to vacuum (the so- 
called standard quantum limit) [2|, |3j, and degenerate op- 
tical parametric oscillators (DOPOs) have been used to 
this aim. A DOPO consists on a nonlinear crystal 
placed inside an optical cavity; when pumped above some 
threshold level with a laser beam at frequency 2lvq, a 



field oscillating at half that frequency {signal frequency) 
is generated. A linearized analysis of the DOPO's quan- 
tum fluctuations reveals that the phase quadrature of the 
signal field can be perfectly squeezed when pump's power 
is tuned to the DOPO's threshold value [J]. Of course, 
ideal perfect squeezing cannot be real in this case, as it 
would entail an infinite number of photons in the gener- 
ated mode, and it can be shown that nonlinear correc- 
tions make this squeezing level become finite 

Apart from fundamental reasons, improving the qual- 
ity of squeezed light is an important task because of its 
applications. For example, in the fields of quantum infor- 
mation with continuous variables [f| (as mixing squeezed 
beams with beam splitters offers the possibility to gener- 
ate multipartite entangled beams Q) and high precision 
measurements (such as beam displacement and pointing 
measurements [|[ or gravitational wave detection 0), ap- 
plications of squeezed light have been theoretically and 
experimentally proved. 

Recently, our group has developed a strategy that 
would lead to the generation of light with high level of 
squeezing at any pump level above threshold in DOPOs 
{non-critically squeezed light). The idea is to allow for 
the existence of several transverse modes at the signal 
frequency in a DOPO possessing some spatial symmetry 
in the transverse plane (e.g., rotational or translational) ; 
under some circumstances, a pattern breaking the corre- 
sponding spatial symmetry can be generated, and quan- 
tum noise can randomly move it along the invariant direc- 
tion (i.e., rotate or translate it in the transverse plane). 
This can be seen as an indefinitcness in the transverse 
position of the generated pattern, which invoking now 
the uncertainty principle, could be accompanied by the 
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perfect determination of its associated momentum (an- 
gular or linear). Hence, we could expect noise reduction 
in the empty pattern coinciding with the momentum of 
the generated pattern. Moreover, as this is related to the 
spontaneous symmetry breaking, the perfect squeezing 
should occur at any pump level above threshold. 

In Ref. [ic| , the translational symmetry breaking was 
considered in broad area, planar DOPOs, where cavity 
solitons (CSs), as well as extended patterns, have been 
predicted to exist [ll[ . It was shown that all the reason- 
ing above is true: The position of the CS diffuses in the 
transverse plane, and the phase quadrature of the pat- 
tern coinciding with its linear momentum (namely its 7r/2 
phase shifted transverse gradient) is perfectly squeezed 
(within the linearized theory) at any pump level above 
threshold. The problem of this model is that it is not too 
close to current experimental setups and even CSs have 
not yet been observed in DOPOs. 

This was one of the reasons why the theory of a 
simpler system was developed in Ref. 12]: A rota- 
tionally symmetric DOPO tuned to the first family of 
transverse modes at the signal frequency [H, [lj|. The 
first transverse-mode family supports Laguerre-Gauss 
(LG) modes L±\ (r) with ±1 orbital angular momentum 
(OAM), where r = (x,y) are the transverse coordinates. 
When pumped with a Gaussian mode with zero OAM, 
two signal photons with opposite OAM are generated in 
the crystal. Alternatively, the simultaneous genera- 
tion of one L + i (r) photon and another L_i (r) one cor- 
responds to the generation of two Hermite-Gauss (HG) 
TEMio photons, the orientation of this TEMio mode {6 
in Fig. 1) given by half the phase difference between the 
subjacent L±\ (r) modes. Again in this case, we proved 
our reasoning given above: Quantum noise is able to ro- 
tate randomly this bright mode, and the phase quadra- 
ture of its angular momentum (its — 7r/2 phase shifted 
angular derivative), which corresponds in this case to 
another HG mode orthogonal to the generated one (we 
will refer to it as the dark mode), is perfectly squeezed 
within the linearized theory and at any pump level above 
threshold [H In addition we showed that this re- 
sult is quite robust against deviations from the perfect 
rotational invariance of the DOPO. 

Moreover, a surprising result was found: although the 
phase quadrature of the dark mode is perfectly squeezed, 
its corresponding amplitude quadrature has not increased 
uncertainty, i.e., the quadratures of this rotating mode do 
not form a canonical pair. We intuitively explained this 
by noting that the excess of noise is transferred not to 
another quadrature, but to the orientation 9 of the pat- 
tern, which actually is the canonical pair of the squeezed 
quadrature; in the present article we prove this statement 
by using the arguments developed in Ref. [l7j . 

Indeed, the study of this rotationally symmetric 
DOPO was carried out under some assumptions whose 
repercussion will be analyzed in the current article. First, 
we supposed that cavity losses at the signal frequency 
were small as compared with that at the pump frequency, 




Figure 1.- Scheme of the 2-transverse-mode DOPO. 



what allowed us to adiabatically eliminate the pump 
variables. Next, the usual linearization of the nonlin- 
ear Langevin equations around their classical stationary 
solution was done. Finally, in order to prove the non- 
critical squeezing we proposed an homodyne detection 
scheme in which the local oscillator was a TEM 10 mode 
matched to the dark mode at any time, i.e., it was follow- 
ing the random rotation of this mode, an ideal situation 
which is not possible in real experiments. 

Now, let us summarize what we are going to show in 
this article, which is divided into four main parts. In the 
first one (Section II) we describe the DOPO tuned to the 
first transverse family (the 2-transverse-mode DOPO) 
and find its associated quantum Langevin equations. In 
Section III, we study the classical emission of the system, 
and show how the rotational symmetry is broken by the 
generated transverse pattern. Section IV is the main sec- 
tion of the article: The quantum properties of the DOPO 
are analyzed by linearizing the Langevin equations. Here 
we will prove the random rotation of the classical pattern, 
and find the squeezing properties of the bright and dark 
modes. In addition, we will show that the quadratures of 
the dark mode do not form a canonical set, but they are 
canonically related to the pattern orientation (where the 
excess of noise goes). We close this section by analyzing 
the repercussion of keeping fixed the orientation of the 
local oscillator, showing that even in this case large levels 
of squeezing can be obtained. In the last part (Section 
V) we show that the approximations considered in the 
previous sections (the adiabatic elimination of the pump 
and the linearization procedure) are not responsible for 
the remarkable squeezing properties of the dark mode, 
as they follow directly from the spontaneous rotational 
symmetry breaking. At the end (Section VI) we will give 
some conclusions. In order to make clearer the physics 
behind the phenomenon, we have left most of the tech- 
nical details to appendices. 
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II. DESCRIPTION OF THE 2-TRANSVERSE 
MODE DOPO 



One can thus define boson operators associated to 
these modes, whose relation with the LG ones is 



A. The field inside the cavity 

One of the simplest configurations of the DOPO's cav- 
ity allowing for the generation of a rotationally asymmet- 
ric pattern is the following [T3|: At some phase-matched 
frequency 2loq it is tuned to a fundamental Gaussian 
transverse mode G(r), while at half that frequency, u>o, 
the first transverse family is resonant. In the first trans- 
verse family two Laguerre-Gauss modes exist, L±\ (r). 
The exact expression of the modes in the waist plane of 
the resonator (where the x^ crystal is placed) is given 
by [3 



G (r) = JlJ-e 

V 7T Wp 

L±i (r) = —j=— je 



(1) 



AM 



where r = r (cos <f>, sin <fi) is the coordinate vector in the 
transverse plane, and Wi is the beam radius at the waist 
plane, which depends on the cavity geometry and the 
resonance frequency of the mode (see Appendix A). 

These transverse modes have and ±1 OAM respec- 
tively, as they are eigenmodes of the OAM operator —id^ 
with those eigenvalues. 

The field inside the cavity can be written then as 

E = iF p A p (r, t) e- 2iW0i + iF s A s (r, t) e~ iuot + H.c, (2) 



with 



(3) 



being n the refractive index of the crystal and L the 
effective cavity length (H.c. stands for "Hermitian con- 
jugate"). The slowly varying envelopes are given by 

A p (r,t) = a (t)G(r) (4a) 

A s (r, t) = a+i (t) L+i (r) + a_i (t) L_i (r) , (4b) 

where the boson operators satisfy the usual equal time 
commutation relations 

[a m (t) , a [ n (£)] = 5 mn ; m, n = 0, ±1. (5) 

Instead of using the LG basis, one can work in the 
most usual TEM m „ (HG) basis. Denoting by H{ (r) and 
(r) a TEMio mode rotated an angle ip with respect 
to the x axis and its orthogonal, these are given by 

Hf (r) = -L [e-**L +1 (r) + e**L. x (r)] (6a) 



= V2|L ±1 (r)|cos (<f>-xP) 
-L[e-*n +1 (r)-e^. 
>/2|£±i(r)|sin (0- i>). 



< (r) = [e-**L + i (r) - e*+L^ (r)] (6b) 



610,^ = — f= (e^a+i + e *^0-i) 



(7a) 
(7b) 



These relations will be important in order to understand 
the properties of the system. 

Finally, note that for any mode m, quadrature can 
be defined as 



X v — e~ itp a -I- e iv n^ 



(8) 



and quadratures with ip = and ip — ir/2 are usually 
called amplitude and phase quadratures, and will be de- 
noted by X and Y, respectively. 



B. Model equations of the system 

Once the structure of the field inside the cavity has 
been described, we can pass to describe its evolution, 
and to this aim a Hamiltonian must be built. It has 
to take into account two processes: the pumping of the 
coherent, Gaussian, resonant laser source at frequency 
2woi and the parametric down conversion of the pump 
photons into signal photons occurring in the x^ crystal. 
In the interaction picture, this Hamiltonian is given by 



H = ih 



Xaoa' +1 a'_ 



H.c. 



(9) 



where £ p and x are proportional to the injected pump 
amplitude and to the second order nonlinear susceptibil- 
ity of the nonlinear crystal respectively. Explicit expres- 
sions of these parameters in terms of physical quantities 
are given in Appendix A. We take £ v as real, i.e., we take 
the phase of the pump laser as the reference phase for all 
the other fields. 

The first part of the Hamiltonian is the usual pump 
Hamiltonian p^ |. while the second part is easily justified 
by energy and OAM conservation: From one 2loq pho- 
ton with zero OAM, two u>o photons are created with 
opposite OAM, or viceversa. Note that this Hamiltonian 
is equivalent to that for the nondegenerate optical para- 
metric oscillator considered for example in Refs. (lq.l20j. 

There exists one more process not taken into account 
in this Hamiltonian: The losses through the cavity mir- 
rors. This is an irreversible process which cannot be de- 
scribed with a Hamiltonian formalism that takes into ac- 
count only the cavity dynamics. However, assuming that 
the outer of the cavity consists of a continuum of modes 
in a vacuum (or coherent) state which weakly interact 
with the internal modes through the partially reflecting 
mirror, it can be incorporated into the master equation 
satisfied by the density operator of the system as shown 
for example in Ref. [21] . Moreover, the master equation 
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can be converted into a set of stochastic (Langevin) dif- 
ferential equations by using a positive P representation 
for the density operator [22j |. In our case, as is shown in 
[20I ] , the equivalent set of Langevin equations is 



mensionless time) 



do = £ p -j p ao - xa+ia-i 

-7 s a+i + x«oaii + s/xao£ (*) 

7 s a_i + x^oaXi + Vxa-oC (t) 
l s at 1 + xaj a+i + \Jx<*o [£ + (*)] ' 



(10) 



ut 

d+i 



where a m and a+ are independent complex amplitudes, 
7 P and 7 S are the decay rates of the cavity at the pump 
and signal frequencies respectively (see Appendix A), and 
the independent complex noises £ and £ + have zero mean 
and nonzero correlations 

(c (to r (k)> = (e + (*i) [e + (*a)] *) = «(*i- * 2 ) . 

(11) 

The equivalence between the master equation of the 
system and these equations must be understood in the 
following way: 



(: / (a m ,a^) :) = (/ (a m ,a+)) 



stochastic 



(12) 



i.e., quantum expected values of normally ordered func- 
tions equal stochastic averages of the same functions 
changing boson operators (& m , a^J by independent com- 
plex stochastic variables (a m ,aj,). Note that these 
stochastic equations are equal in either the Ito or 
Stratonovich forms (see Appendix B), and hence we in- 
terpret them as Stratonovich equations which allow us to 
apply ordinary calculus. 

In order to better visualize the free parameters of the 
model, let us define new rescaled time and amplitudes 
through 

t = 7s*) 0+1 = — a±i, 0a = — ato, (13) 

VTpTs 7s 

and similar expressions for the conjugate fields aj^ and 
Oq , Note that as time is measured in units of 7s -1 , the 
noises must be rescaled also as 



C±i (r) = — £ ±x (t) , 
\/7s 



(14) 



and similarly for £± lt if we want to preserve their sta- 
tistical properties (jlip in terms of the adimensional time 
r. 

Finally, these changes make Langevin equations (I10p 
read (of course, derivatives are respect with the new di- 



/?o = 


K ((J — 


- fa - 0+10-1) 




4 + = 


K (<7 — 


- 0X - 0Xi0U 


) 


P+i = 


-0+1 


+ 0a0-i +9V 




0\i = 


-0Xl 


+ 0X0-1+9X 


/^C + (*) 


0-i = 


-0-1 


+ O 0+ 1 + <?V 


^C* (*) 


PU = 


-0\ 


+ (3X0+1 + 9\ 


/^[C + (<)]* 



(15) 



which have only the following 3 dimensionless parame- 
ters: 

7 P £ P x 
k = — , a = — — 

7s 



and q = 



X 



7 P 7s 



III. CLASSICAL EMISSION: 
SYMMETRY BREAKING 



VTpTs 



(16) 



ROTATIONAL 



Before analyzing the quantum properties of the sys- 
tem, let us examine the emission of this DOPO as pre- 
dicted by classical optics. It is possible to retrieve the 
classical equations of the 2-transverse-mode DOPO from 
the quantum Langevin ones (fH)]) by setting all the noises 
to zero and making /?+ = (to refers to any mode of 
the problem) . We are interested in finding the stationary 
emission, and then set to zero the time derivatives of the 
amplitudes, arriving to 

Q = a- /3+i/3_i (17) 

+1 = O 0*_ X 

0^=0 O 0* +V 

Before solving these equations, let us stress that these 
are invariant under the following transformation 







±1 



±ii/> 



0+1 



(18) 



and hence the phase difference between opposite OAM 
modes is classically undetermined. 

By decomposing into modulus and phase the ampli- 
tudes m , it is straightforward to find the solutions of 
the system (IT7)) . If in addition a linear stability analysis 
is made, it is easy to find the following result: 



• For a < 1 the only stable solution is 

0o = @ 
0+1 = 0, 



(19) 



and hence the signal modes are off if the pump £ p 
doesn't exceed a threshold value Eth = 7p7s/x- 

• On the other hand, if a > 1, the signal modes are 
switched on, and the only stable solution in this 
case is 



00 = 1 
0+i = pe Ttd with p : 



(20) 



where, as stated above, 9 is arbitrary. 
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If we substitute the last result for the signal modes in 
the expression of the corresponding slowly varying enve- 
lope (|4b]l. we find the generated pattern to be 



A s (r) = p [e- w L +1 (r) + e+ ie L_ 1 (r)] = yfipH% (r) , 

which is a TEM 10 mode forming an angle 9 with respect 
to the x axis (O as shown in Fig. 1. Due to the arbi- 
trariness of 9, the pattern can arise with any orientation; 
this reflects the rotational invariance of the system, wich 
in turn is broken after the TEMio mode generation. 

In the following, we will call Hf (r) the bright mode (as 
it is classically excited) and its orthogonal mode H® x (r) 
the dark mode (as it is classically empty of photons), 
and will define the collective indices b = (10, 9) and d = 
(01,9) to simplify the notation. 



IV. QUANTUM PROPERTIES: PATTERN DIF- 
FUSION AND NON-CRITICAL SQUEEZING 

A. The linearized Langevin equations 

We are going to discuss the quantum properties of 
the down converted field by inspection of the quantum 
Langevin equations (|15[) in the limit 7 P 2> 7 S , i.e., k» 1, 
where the pump variables can be adiabatically elimi- 
nated. We will show later (Section V) that all the impor- 
tant properties found in this limit are valid in general. 

One common way of treating the adiabatic elimina- 
tion of the pump consists on setting (3q — 0q — in the 
first two equations of ([15)) . However, this method has a 
problem: The initial Langevin equations are equal in ei- 
ther Ito or Stratonovich forms, but the new Langevin 
equations obtained for the signal modes doesn't have 
this property (see Appendix B). Hence it is fair to ask 
within which interpretation (Ito or Stratonovich) is that 
procedure correct, if it is correct at all. As proved 
in (23j by performing the adiabatic elimination in the 
Fokker-Planck equation (where there are no problems 
of interpretation), the usual method for eliminating the 
pump is correct within Ito's interpretation, and hence 
the Stratonovich form of the stochastic equations satis- 
fied by the signal modes only has an unexpected extra 
term proportional to g 2 (see Appendix B) and reads: 



= 

0tx = 
$-1 = 

fcx = 

with 



-(1- 3 2 /4)/? +1 +/3 /3± 1 +. 9 V^C(t) (22) 

- (1 - , 9 2 /4) p+i + tiP-i + 9\[tiQ + (r) 

- (1 - g 2 /4) /3_ a + poPti + 9VWoC (r) 

- (1 - 5 2 /4) (3± x + + gJ(3+ [(+ (r)] * , 



00 = a - /3 + i/3 + i, 
/? + = a - PtiP+i- 



(23) 



In order to find analytic predictions from these equa- 
tions, we assume well above threshold operating condi- 
tions, where classical emission dominates, thus allowing 
for a linearization procedure. The correct way to perform 
this was described in [f| [24| , and relies on the smallness 
of the g parameter (see Appendix A). In particular, mak- 
ing a perturbative expansion of the amplitudes in terms 
of this parameter, it is possible to show that the classical 
equations are recovered as the g° order of the expansion 
in Eqs. (|22p . while the linear quantum correction ap- 
pears in the g 1 order (note that the noise term is already 
of order g 1 ). 

Hence, the usual linearization procedure begins by 
writing the amplitudes as (3 m — /3 m + Sf3 m and /3+ = 
/3* t + 5/3+ , and treat the fluctuations as order g pertur- 
bations. However, as stated in the Introduction, in the 
system we are dealing with we expect that quantum noise 
rotates the generated TEM10 mode, and hence fluctua- 
tions of the fields in an arbitrary direction of phase space 
could not be small (i.e., order g). Nevertheless, Eqs. 
can be linearized if the amplitudes are expanded as 



/3± 1 (r) = [p + 6 ±1 (r)]e = F < °W 
Pti (r) = [p + b+ 1 (r)] e^\ 



(24) 



because as we will prove 9 (r) carries the larger part of 
the fluctuations, while the &'s and 9 remain as order g 
quantities. In addition, expanding the fields in this way 
allows us to track the evolution of the classical pattern's 
orientation, as we take 9 as an explicit quantum variable. 
Then, writing Eqs. (|22|) up to order g, we arrive to the 
following linear system (arriving to this expression is not 
as straightforward as it might seem, see Appendix C) 



2ipw <9 + b =Ch + gC (r) , 



with 







f C(r) \ 


b-i 


,C(r) = 


C + (r) 
C* (r) 






\K + (r)]*J 



(25) 



(26) 



and where £ is a real, symmetric matrix given by 



C 



with the following eigensystem 



a 





a- 1 


-1 





a 


-1 


a — 


- 1 


-1 


a 





-1 


a - 1 





a 



(27) 



A,, 


= 0, 


w 


Ai 


= -2, 


Wl 


A 2 


= -2(ct- 


1) , w 2 


A3 


= -2a, 


W3 



§001(1,-1,-1,1) 



\ col (1,1, -1,-1) 

2 = I COI (1,1, 1,1) 

3 = 1 col (1,-1, 1,-1) 



(28) 



Defining the projections c m (r) = w TO • b (r), and pro- 
jecting the linear system (|2~5j) onto these eigenmodes, we 
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find the following set of decoupled linear equations (co is 
set to zero, as otherwise it would just entail a redefinition 
of 6) 

0=f p V° (r) (29a) 

ci = -2ci + ig m (r) (29b) 
c 2 = -2 (<r - 1) c 2 + ,g?72 (r) (29c) 
c 3 = -2c7c 3 + 3773 (r) , (29d) 

where the following real noises have been defined 

, 7o (r) = zw • C (r) = Im { (+ (r) - C (r) } (30) 
m (r) = -t Wl • C (r) = Im {C+ (r) + C (r)} 
^2 (r)=w 2 .C(r)=Rc{C(r)+C + (r)} 
m (t)=w 3 -C(t)=Rc{C(t)-C + (t)}, 
which have zero mean and nonzero correlations 

(Vm ( T i) ih. (r 2 )) = S mn S (n - r 2 ) . (31) 

Note finally that the solutions for the projections cj 
will be of order g (see Appendix D), and hence so will 
be the b's (note that this is not the case for 8, whose 
initial value is completely arbitrary, although its varia- 
tion is indeed of order g) . This is consistent with the 
initial assumptions about the orders in g of the involved 
quantities. 

B. Quantum diffusion of the classical pattern 

The first quantum effect that we are going to show con- 
cerns the orientation of the classically generated mode. 
Equation (|29a[) defines a Wiener process for 8, thus 
showing that the orientation of this bright mode diffuses 
with time ruled by quantum noise. 

How fast this diffusion is can be measured by evaluat- 
ing the variance of 8. Using the statistical properties of 
noise pip , it is straightforward to obtain the following 
result 

V e (T) = (58 2 (T)) = DT, (32) 
where D = d/ (a — 1) with 

d = 5 2 /4 = X 2 /47 P 7s- (33) 

In P"2"| the notation SA = A — (A) has been used. 

Eq. (J2U) shows that the delocalization of the pattern's 
orientation increases as time passes by, though for typical 
system parameters (see Appendix A) one finds d ~ 4 • 
1CP 13 , and hence the rotation of the pattern will be fast 
only when working terribly close to threshold [25| . 

C. Squeezing properties of the 2-transverse-mode 
DOPO 

Definition and criterium for squeezing. In the In- 
troduction we defined squeezed light as that having some 



quadrature fluctuating below the vacuum level. Hence, 
in order to find out whether a light beam is in a squeezed 
state or not, one has to measure its quadrature fluctua- 
tions and then compare them with the value set by the 
vacuum state. The question then is: What quantity ac- 
counts for quadrature fluctuations of a light beam? Of 
course, in the case of a single mode of light this quan- 
tity can be directly the uncertainty of the quadrature. 
However, outside the cavity there exists a continuum of 
modes, and the quantity accounting for these fluctuations 
has to be adapted to what can be most easily observed 
in an experiment. 

As was first shown in [J], the quantity accounting for 
fluctuations of quadrature X^ outside the cavity (m 
refers to any signal spatial mode of our 2-transverse-mode 
DOPO) is 

V out (u>;Xv) =1 + S*(u), (34) 

with 

S ™ M = 4 df (: 8X% (r) 6X% (r + f) :) 

9 J — OO 

(35) 

where the factor g 2 appears after including our rescaled 
variables (fl~3|) . 

We will call V ont the noise spectrum and 5^ (u>) the 
squeezing spectrum. Frequency u> is usually called noise 
frequency, and it must not be confused with the optical 
frequency. In fact, noise frequency uj has contributions of 
every pair of modes lying in opposite sidebands around 
the optical frequency u)q + u>, where luq is the carrier fre- 
quency of the detected beam (don't miss that in our case 
this frequency is measured in units of 7 S ). 

V out can be measured via a balanced homodyne detec- 
tion experiment: A coherent, intense field (local oscillator 
field) prepared in mode to, and with a phase (p is mixed 
with the beam exiting the DOPO in a 50/50 beam split- 
ter; it is easy to show that the operator associated to the 
intensity difference between the two output ports of the 
beam splitter is proportional to the quadrature X^ ont of 
the beam exiting the DOPO. Hence, the normalized cor- 
relation spectrum of the intensity difference (which can 
be measured with a simple spectrum analyzer) coincides 
with expression (|34"|) . Factor 2 in ([35]) comes from the 
relation between the intracavity and the output modes 
when the input modes (the pump laser and the outer 
modes) are coherent; this factor is 2^ s when time (fre- 
quency) is not normalized to 7 s _1 (7s)- 

In the above expression for the squeezing spectrum, 
it is assumed that the emission of the DOPO reaches a 
stationary state (i.e., the correlation function of SX^ at 
two different times r and r + f depends only on the time 
difference |f|) and that the observation time, say T, is 
large as compared with the coherence time of SX^. In 
the 2-transverse-mode DOPO there exists an undamped 
quantity excited by quantum noise, 9, and hence, any 
quantity with explicit dependence on it will not arrive to 



G 



a stationary state. Thus, under some circumstances we 
will need the expression of the squeezing spectrum mea- 
sured when the field is not stationary, which was derived 
in 26]. In that reference it is shown that when the mea- 
surement extends during a time interval T (normalized 
to 7,7 1 in our expressions) the squeezing spectrum reads 



si M = 



9 2 TJ 



drdr' (:5XI(t)5XI(t>): 



x cos [uj (t - r')] , (36) 



which obviously reduces to (|35[) if the emission is station- 
ary. 

For the vacuum state V out — 1, hence if 

V° ut (lj: X^j < 1 for some frequency u>, we can say that 

the detected light is in a squeezed state for mode m at 
noise frequency u>. 

Finally, let us stress that we will evaluate quantum 
correlations as those involved in (|36p via stochastic cor- 
relations, as allowed by property ([12"]) . Hence, in the 
following we will change operators by their equivalent 
stochastic variables within the positive P representation 
in all the definitions above, and their associated normally 
ordered correlations by the corresponding stochastic cor- 
relations. 

Independent quadratures and non critical 
squeezing. A first step towards analyzing the squeezing 
properties of the field within the framework presented 
above, is identifying a set of independent quadratures. 
This is easy in our case, as the eigenmodes w m give us a 
set of quadratures with well defined squeezing properties. 
In particular, from ([24]) . ([7]), and ([8]), it is easy to find 
the following relations 



X h (t) = 2V2p + V2c 2 (r) (37) 

Y b (r) = -iV2c 3 (r) 
X d (r) = iV2c (r) 

Y d (r) = y/2d (r) , 

where (X\,, It,) and (Ad, Y d ) are the amplitude and phase 
quadratures of the bright and dark modes, Hf Q (r) and 
#oi ( r )j respectively. 

The evolution of these quadratures can be found from 
the equations satisfied by the projections (|29|) . which are 
solved in Appendix D (remember that cq = 0). In partic- 
ular, there we show that after some time these projections 
reach a stationary state, and hence we can evaluate the 
noise spectrum of the quadratures by using the station- 
ary expression for the squeezing spectrum (|35|) . Using 
the results found in Appendix D for the correlation spec- 
trum of the projections (|D7[) . it is straightforward to find 



the following results 



V mt (cu:X k 



1 



V out (w;y b ) = l- 
V° ut (w;X d ) = l 
V mt (w;Y d ) = l- 



{a-iy+u 2 /4, 
1 



a 2 + uj 2 /4 
1 

1 +w 2 /4' 



(38a) 

(38b) 
(38c) 
(38d) 



We see that the quadratures of the bright mode have 
the same behavior as those of the single mode DOPO: 
the phase quadrature Yb is perfectly squeezed (V out = 0) 
at zero noise frequency {uj = 0) only at the bifurcation 
(a=l). 

On the other hand, the dark mode has perfect squeez- 
ing in its phase quadrature Yd at zero noise frequency. 
What is interesting is that this result is independent of 
the distance from threshold, and thus, it is a non-critical 
phenomenon. 

This result was first shown in [12j, and we can under- 
stand it by following the reasoning given in the Intro- 
duction: As the orientation 9 of the classically excited 
pattern is undetermined in the long-time limit, its OAM 
must be fully determined at low noise frequencies. On 
the other hand, the OAM of the bright mode Hf Q (r) is 
nothing but its ir/2 phase shifted orthogonal HG mode, 
i.e., —id<f,Hf Q (r) = iH^ (r), which used as a local oscil- 
lator in an homodyne detection experiment would lead 
to the observation of the Yd quadrature fluctuations. 

On canonical pairs and noise transfer. Although 
all the results we have discussed up to now were expected 
to occur after the arguments given in the Introduction, 
a strange unexpected result has appeared in ([38)) : The 
quadratures of the dark mode seem to violate the uncer- 
tainty principle as V out {uj = 0; X d )-V out (uj = 0; Y d ) = 0. 

Moreover, the noise spectrum of an arbitrary quadra- 
ture of the dark mode, which can be written in terms of 
its amplitude and phase quadratures as X^ — X d cos ip + 
Y d sin ip, is 



Y out 



1 



sin 2 ip 



(39) 



what shows that all the quadratures of the dark mode 
(except its amplitude quadrature) are squeezed. Thus 
two orthogonal quadratures cannot form a canonical pair 
as they satisfy the relation 



V out {u;X%)-V out (w;XX +7r/2 ) 



< 1, 



(40) 



in clear violation of the uncertainty principle. 

The natural question now is: Where does the excess 
of noise go if it is not transferred from one quadrature 
to its orthogonal one? The answer is that it goes to the 
pattern orientation, which is actually fully undetermined 
in the long term as we showed above This section 

is devoted to prove this statement. 
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In particular we will prove that two orthogonal quadra- 
tures of the dark mode do not form a canonical pair, 
while the orientation 9 is the canonical pair of all the 
squeezed quadratures. One way to prove this would be 
to evaluate the commutator between the quantum op- 
erators involved. However, 9 is half the phase difference 
between the opposite OAM modes /3±i, whose associated 
operator has a very difficult expression [2?| . making the 
calculation of the needed commutators quite hard. Nev- 
ertheless, in [13] we developed a simple proof based on 
classical field methods, and this is the one we present 
here. 

A usual approach one uses to move from classical to 
quantum optics is to change the classical normal variables 
for each mode of the field, (3 m and [3 m , by boson operators 
a m and a m , satisfying commutation relations 



'mi &m] 



i {Pirn Pm} > 



(41) 



where {/, h} denotes the Poisson bracket between two 
functions / (/3 m , and h ((3 m , /3£J defined as 

u ' 1 i^dPmdfa dp* m dp m v ; 

As an example, the Poisson bracket of two monomode 
orthogonal quadratures X m = (3 m + f3 m and Y m 
i (/3 m — /?£J is found to be 



(43) 



In general two phase space functions /(/3 m ,/3^J and 
h {(3 m , P m ) are said to form a canonical pair if their Pois- 
son bracket is of the form 



{f,h} = C, 



(44) 



with C a real number. 

In our case, the functions we are interested in are 
the classical counterparts of the dark mode quadratures, 
which using (|7b|) with -0 = 9 are written as 

XI = -L [e-* (e*/? +1 - e- 1 ^)] + c.c, (45) 



with the orientation given by 



2i 



In 



^PUPVi^PUP- 



(46) 



Now using the definition of the Poisson brackets with 
m = ±1 and after some algebra, it is possible to show 
that 



and 



\P-i\-\P+i\ 

2 |/3 + i I l/3-i I ' 
\P-i\) (e^ \P-i\ 



(47) 



-tif 



1/5 1 1) 



4V2|/3 +1 | 3/2 |/3_ 1 | 



,3/2 



(48) 



With \f3 m \ = \/P* m P m - 

On the other hand, in the two-transverse-mode DOPO, 
the number of photons with opposite OAM is sensibly 
equal, i.e., ~ \P-i\'i hence, the dominant term of 

the previous brackets will be that with |/3+i| = \P-i\ = p, 
and thus 



XT 



-7T/2 







and 



siny 

Tip' 



(49) 



(50) 



Comparing these with (|4"4"|) we find that this is a con- 
firmation of what we expected: two orthogonal dark 



quadratures Xf and X 



V+n/2 



do not form a canonical 



pair (moreover, they commute), while X^ and 9 do. This 
can be seen as an indirect proof of the same conclusion 
for the corresponding operators. 

Homodyne detection with a fixed local oscilla- 
tor. So far we have considered the situation in which 
one is able to detect independently the bright and dark 
modes. However, as shown by Eq. (|29a[) . these modes 
are rotating randomly, what means that a local oscilla- 
tor field following that random rotation should be used in 
order to detect them separately. This might be a really 
complicated, if not impossible, task, so we analyze now 
the more realistic situation in which the local oscillator 
is matched to the orthogonal orientation of the emerging 
pattern only at some initial time, remaining then with 
the same orientation during the observation time. We 
will show that even in this case, and as the rotation of 
the modes is quite slow (|3"2j) . large levels of noise reduc- 
tion can be obtained. 

Without loss of generality, we suppose that the bright 
mode emerges from the resonator at some initial time 
t = oriented within the x axis, i.e., 9 (0) = 0. Hence, 
by using a TEMoi local oscillator with a phase <p, the 
quadrature X^ of a fixed H i (r) mode (the initially dark 
mode) will be measured. In terms of the Gauss-Laguerre 
modes, the amplitude /3 i of the this mode is given by 
(|7b[) with ip = 0. Then, by using the expansion of 
the amplitudes (3 m as functions of the fluctuations b m 
and the orientation 9, quadrature Xq X of this mode can 
be rewritten as 

X& = 2V2p cos ip sin 9 + V2c 2 cos(ps'm9 (51) 
+ V2ci sin tp cos 9 — i\picz sin ip sin 9. 

Hence, the 2-time correlation function of X^ yields (note 
that (f2"9"]) shows that c\ , Ox , C3 , sin 9 and cos 9 are uncor- 
rected) 

(X»' 1 (t 1 )X ¥: 1 (t 2 )) = 2cos 2 pS (7-1,7*) [4p 2 + C 2 (7i,T 2 )] 

(52) 

+ 2 sin 2 p [d (n, 75) C (ti,t 2 ) - C 3 (71,7-2) S (7-1,7-2)] , 



8 



with 



S (ri,r 2 ) 

C(t u t 2 ) 

C m (7-1,75) 



(sin^ (n) sin^ (72)) 
(cos 9 (ti) cos (9 (72)) 
(c m (n) c m (r 2 )) . 



(53) 



As shown in Appendices D and E, the latter correlation 
functions can be evaluated by using the linear evolution 
equations of the projections c m and 9 -see Eqs. (|D6I) . 
(|E4j) and (|E5| -: then from (|52|) the squeezing spectrum of 



Xq^ can be found by using the general expression (f3"6"|) , as 
in this case the stationary expression ([33)) cannot be used 
because S* (71,72) and C (71,72) don't reach a stationary 
state. 

Although the general expression for Sq X (lj) is too 
lengthy to be written here, a more compact approximated 
expression can be found in the limit of small d, leading 
to the following expression for the noise spectrum 

V mt (w; X&) = (w) cos 2 cp+£%{ 2 (w) sin 2 (54) 



with 

Sgi = -^(1-sincwT)- 



4dT 6(cr-l) 2 



w 2 (ct-1) 4(cr-l) 2 +u 2 



and 



Q 7T/2 _ 

°oi — 



(55) 



■ -2w 2 



8cff [2 (cr 2 + 1) + to 2 ] 



T{A + lo 2 ) 2 4 + w 2 ' ( ( x-l)(4 + ^ 2 )(4a 2 +w 2 )' 

(56) 

where sine x — sin (x) jx. In the following we will fix a to 
\/2 (pump power twice above threshold), as the results 
are almost independent of its value as far as it is far 
enough from threshold. Hence, the free parameters will 
be the detection parameters T, uj and tp, and the diffusion 
d which depends on the system parameters. 

In this section, the results for V out are presented 
in dB units, defined through the relation V out [dB] = 
101ogy out (hence, e.g., —10 dB and —00 dB correspond 
to 90% of noise reduction (V out = 0.1) and complete 
noise reduction (V out — 0) respectively.) 

From expression (|54p we see that the maximum level 
of squeezing is obtained at u> = (see also Fig. 3a) and 
when the phase of the local oscillator is tuned exactly 
to 7r/2. In Fig. 2 we show the noise spectrum for 
these parameters as a function of the detection time T 
for 3 different values of d. We see that in all cases there 
exist an optimum detection time for which squeezing is 
maximum. Minimizing Eq. (|54f> with u> = and tp = tt/2 
with respect to T, it is straightforward to find that this 
optimum detection time is given by 



T, 



opt 



I a 2 (<7-l) 
d(a 2 + l)' 



with an associated noise spectrum V^° ut 



v opt 



in the inset of Fig. 2 as function of d). 



(57) 



l/^opt (shown 



These results show that large levels of noise reduction 
are obtained for the phase quadrature of the fixed TEM i 
mode, even for values of the diffusion parameter d as large 
as 10~ 6 (remember that 10~ 13 is a more realistic value). 
However, in real experiments it is not possible to ensure 
that tp = 90 with an uncertainty below approximately 

i.5° nm, and hence we proceed now to investigate the 

level of noise reduction predicted by l|54p when the local 
oscillator phase is different from tp = 90 . 

Of course, when ip ^ 90° the noise frequency with 
maximum squeezing is no longer lo = 0, as in this case 
the infinite fluctuations of at zero noise frequency, 
due to the rotation noise, enter the noise spectrum (see 
Fig. 3a). By minimizing y ont with respect to u and 
T for different values of tp and d, it is possible to show 
that the optimum value of the detection time is almost 
independent of tp for small deviations of this from 90°, 
and hence it is still given to a good approximation by 
(|57p , though in this case this minimum is less pronounced 
than in the tp = 90° case shown in Fig. 2 (i.e., the curve 
is almost horizontal around T opt ). On the other hand, 
the optimum noise frequency w opt is independent of d 
and depends on tp as shown in the inset of Fig. 3b. 

As for the squeezing level, in Fig. 3b we show the 
noise spectrum evaluated at T opt and w op t as a function 
of tp for 2 different values of the diffusion d. Together 
with these curves, we have plotted the noise spectrum 
of the single-mode DOPO [28| evaluated for its optimum 
parameters (in this case it is optimized respect to a and 
to) as a function of tp. We see that the noise reduction is 
independent of d as tp is taken apart from 90 . On the 
other hand, the squeezing level is similar to that of the 
single-mode DOPO, as the maximum difference between 



V out (m = ; Xlf°° ) [dB] 




5 10 

T [dO" 6 ] 

Figure 2.- Zero-frequency noise spectrum of the phase 
quadrature corresponding to a fixed TEM01 mode as a 
function of the detection time T. Three different values of d 



are considered (10 



10" 



and 10 from top to 



bottom). The inset shows also this spectrum at zero noise 
frequency, but evaluated at the optimum detection time T opt 
and as a function of d (note that the d axis is in logarithmic 
scale) . As mentioned in the text a — v2- 



9 



them are 1.5 dB (a factor 1.4 in the noise spectrum) in 
favor of the single-mode DOPO, with the advantage that 
in the 2-transverse-mode DOPO this level is independent 
of the distance from threshold. 

Therefore we see that the phenomenon of non-critical 
squeezing through spontaneous rotational symmetry 
breaking could be observed in the 2-transverse-mode 
DOPO without the need of following the random rotation 
of the generated pattern, which makes its experimental 
realization feasible with current available technology. 



25 



(a) 



f -- 





v\ s 9 = 89 5 ° 


-0.5 




VOs/ 0.5 




N^25 






-50 
v opt 


[dB] 



CO 




-10 



-30 



-50 



90° 



Figure 3.- (a) Noise spectrum of the fixed TEMoi mode as a 
function of the noise frequency u, and for four different 
values of the local oscillator phase (cp — 88 °, 89 °, 89.5 °, 
and 90° from top to bottom curves). It can be appreciated 
how the infinite fluctuations of Soi at zero noise frequency 

enter in the spectrum for any ip 7^ 90. The rest of 
parameters are a = v2, d = 10 -10 , and T = T opt , but the 
same behavior appears for any election of the parameters, 
(b) Noise spectrum of the fixed TEM01 mode evaluated for 
the optimum parameters ui op t and T op t as a function of ip 
(d = 10 -13 for the blue-solid curve and 10 -6 for the 
red-dashed one, having both a = y2). In addition, it is 
plotted the analogous curve for the single-mode DOPO 
(gold, dashed-dotted curve). The inset shows the 
dependence of the optimum frequency w pt with the phase of 
the local oscillator ip. 



V. BEYOND THE CONSIDERED APPROXIMA- 
TIONS 

A. Beyond the adiabatic elimination of the pump 

The first assumption made in the search for the quan- 
tum properties of the system was that 7 P 3> 7 S , a limit 
that allowed the adiabatic elimination of the pump field. 
Now we are going to show analytically that the phe- 
nomenon of squeezing induced by spontaneous rotational 
symmetry breaking is still present even without this as- 
sumption, but still working in the domain where lin- 
earization is correct. 

The way to show this is quite simple; starting from 
the complete equations (fl~5|). we expand the amplitudes 
(3 m around the classical stationary solution (f20|) as we 
made in (|24p , but adding now a similar expression for 
the pump amplitudes: /3q = l + &o and @q = l + b$ (note 
that for the pump modes the 6's are directly small as the 
phase of this mode is locked to that of the injection £ p ). 
Then, linearizing these equations for the fluctuations and 
noises, we arrive to a linear system formally equal to (|25[) . 
but with 



(58) 













f s ) 


b+i 

Ki 


,c = 


CM 

C + (r) 


b-i 




C* (r) 


\b±J 




V[C + (r)]7 



and a linear matrix 



C = 



f-K 





-p 





-p 


°\ 





— K 





-p 





-p 


p 





-1 








1 





p 





-1 


1 





p 








1 


-1 





V 


p 


1 








-l) 



(59) 



Although in this case C is not hermitian, it can be 
checked that it possess a biorthonormal basis [2{|, and 
in particular, the following two vectors are present in 
its eigensystem: Wg = ^ col (0, 0, 1,-1,-1, 1) and w[ = 
i col (0,0, 1,1,— 1,-1), with corresponding eigenvalues 
A = and = —2. These eigenvectors have null pro- 
jection onto the pump subspace, and coincide with w 
and Wi in what concerns to the signal subspace (|28[) . 
Hence, all the properties derived from these vectors are 
still present without any change. In particular, as they 
account for the diffusion of 9 and the squeezing properties 
of the dark mode, we can conclude that these properties 
are still present when working out of the limit 7 P 3> 7 S . 



B. Numerical simulation of the nonlinear equa- 
tions 

In this section we will show that the diffusion of the 
orientation and the associated non-critical squeezing of 
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the dark mode, which have been found by linearizing 
the Langevin equations, are also present when we con- 
sider the full nonlinear problem. To do so, we will solve 
numerically the complete stochastic equations (1151) us- 
ing the semi-implicit algorithm developed by Drummond 
and Mortimer in [30j ] . 

The details of the numerical simulation are explained 
in Appendix F. Here we just want to point out that the 
important parameters of the simulation are the step size 
At used to arrive from r = to the final integration 
time r en( j, and the number of stochastic trajectories, say 
E, which are used to evaluate stochastic averages. The 
initial conditions (3 m (0) are not important as the results 
in the stationary limit are independent of them. The 
system parameters which have been chosen for the simu- 
lation are a = \/2, n — 1 (to show also numerically that 
the adiabatic elimination has nothing to do with the phe- 
nomenon), and g = 1CP 3 . We haven't chosen a smaller 
value for g (like 10 -6 as followed from the physical pa- 
rameters considered in Appendix A) because such a small 
number can make the simulation fail; nevertheless all the 
results we are going to show should be independent of g 
and cr, and we have also tested that the same results are 
obtained for other values of these. 

It is also important to note that we have defined a 
general quadrature of the rotating dark mode as (directly 
from ([7b| with (3 = 6) 



-it) 



0-0] 



(60) 



with 9 defined within the positive P representation 
through 



„2i0 



+i 



(61) 



+i| 



Now let us show the results evaluated for the following 



simulation parameters: At = 3-10 



-3 



Tend 



= 30 and 



£ = 10 7 . This simulation has been compared with other 
ones having different values of these parameters to ensure 
convergence. 

In Fig. 4a we show the variance of 9 as a function of 
time. The result has been normalized to D, so that the 
linear result (j3"2")l predicts a straight line forming 45° with 
respect to the time axis. It can be appreciated that this 
is indeed what shows the simulation. 

In Fig. 4b, we show the numerical results for the noise 
spectrum associated to the quadratures of the dark mode. 
Only times above r = 10 have been considered in the 
correlation function to ensure being working in the sta- 
tionary limit. Again, the results shown in Fig. 4b are in 
perfect agreement with the linear predictions (|38cp and 
pd) . 

These results show that the phenomenon of non- 
critically squeezed light via spontaneous rotational sym- 
metry breaking is not a product of the linearization. 
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Figure 4.- (a) Evolution of the variance of the orientation 6 
given by the numerical simulation. It has been divided by 
the slope D predicted by the linearized theory (|32p . so the 
straight line obtained is in perfect agreement with this linear 
result. Quantitatively, a linear regression analysis shows 
that the slope obtained by the data is 1.00108, with 
standard error 9 • 10~ 5 . (b) Noise spectrum of the amplitude 
(red-dashed curve) and phase (blue-solid curve) quadratures 
of the dark mode as obtained by the numerical simulation. 
The results are in perfect agreement with the ones predicted 
by the linearized theory (|38[) . A nonlinear regression 
analysis of the numerical data respect to the function 
a (lo/2) 2 /[b + a (uj/2) 2 } , shows that (a, b) = (0.9446, 0.9591) 
are the best fit parameters having both standard errors 
4- 10~ 4 , which is in good agreement with the values 
(a, b) = (1, 1) predicted by (J38J). 



VI. CONCLUSIONS 

In conclusion, we have developed in detail the quan- 
tum theory of the 2-transverse-mode DOPO introduced 
in [l2| . We have studied some important features not 
considered in that reference: 

• When pumping with a gaussian mode a rotationally 
symmetric DOPO properly tuned to the first family 
of transverse modes at the signal frequency, classi- 
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cal emission takes place in a TEMio mode with 
an arbitrary orientation in the transverse plane. 
Hence, once the threshold is crossed, the rotational 
symmetry is spontaneously broken, and we can 
talk about a bright mode (the generated one which 
breaks the symmetry) and a dark mode (the mode 
orthogonal to the generated one). 

The symmetry breaking reflects in the appearance 
of a Goldstone mode in the matrix governing the 
linear evolution of the fluctuations above thresh- 
old. The null eigenvalue of this mode allows quan- 
tum noise to change the orientation of the bright 
mode randomly. Though continuously increasing 
with time, this rotation of the classically excited 
pattern is quite slow when working above thresh- 
old. 

As for the squeezing properties, it has been proved 
that the bright mode has the same behavior as the 
single-mode DOPO, i.e., perfect squeezing appears 
only at threshold (within the linearized theory) and 
degrades fast as pump is moved apart from this 
level. On the other hand, accompanying the Gold- 
stone mode it appears another mode whose associ- 
ated eigenvalue takes the minimum possible value 
(—2). These modes are responsible of the remark- 
able properties of the dark mode: Its phase quadra- 
ture is perfectly squeezed at any pump level, while 
its amplitude quadrature carries only with vacuum 
fluctuations [12[ (in apparent violation of the un- 
certainty principle). A simple explanation of this 
phenomenon in terms of " angle - angular momen- 
tum" uncertainty relation appears once one notices 
that the dark mode coincides (up to a tt/2 phase) 
with the OAM of the generated pattern. 

We have proved that the apparent violation of the 
uncertainty principle is just that, apparent, as the 
conjugate pair of the squeezed quadrature is not an- 
other quadrature but the orientation of the bright 
mode, which in fact is completely undetermined in 
the long term. 

Next we have pointed out that in order to mea- 
sure the quantum properties of the dark mode, one 
has to use a TEMio local oscillator that is per- 
fectly matched to the orientation of this mode at 
any time. However, the mode is rotating randomly, 
which seems to make impossible the perfect match- 
ing. For this reason we have studied the situation 
in which the local oscillator is matched to the dark 
mode's orientation only at the initial time, remain- 
ing fixed during the detection time. We have shown 
that arbitrarily large levels of noise reduction can 
be obtained even in this case if the phase of the local 
oscillator is exactly tt/2. We then considered phase 
deviations up to 2° (1.5° seems to be the current ex- 
perimental limit 0, HJ), comparing the results with 
that predicted for the single-mode DOPO; similar 



levels are obtained for both, with the advantage 
that in the 2-transverse-mode DOPO this level is 
independent of the distance from threshold, and 
hence, non-critical. 

• In the last part of the article we have shown that 
the assumptions made in order to analytically solve 
the problem are not the responsible for the quan- 
tum properties of the dark mode. In particular, 
we have shown that when the pump is not adia- 
batically eliminated, the Goldstone mode and its 
companion with the lowest possible eigenvalue re- 
main unchanged in the matrix governing the lin- 
ear evolution of the full problem, and hence all the 
properties derived from them are still present. Fi- 
nally, we have used numerical simulations to show 
that the predictions of the linearized equations are 
in perfect agreement with those of the full nonlin- 
ear equations, and hence the perfect, non-critical 
squeezing of the dark mode is not a by-product of 
the linearization. 

We believe that the analyses presented in this paper 
(together with the fact that the properties of the dark 
mode are not too sensitive to imperfections in the ro- 
tational symmetry of the DOPO [12], and that the phe- 
nomenon is present in other kinds of nonlinear resonators 
[3l| ) show that that the phenomenon of non-critical 
squeezing induced by spontaneous rotational symmetry 
breaking is a robust phenomenon. 
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APPENDIX A: MODEL AND PHYSICAL PA- 
RAMETERS 

During the article some parameters have been used 
to model the DOPO. Here we want to give explicit ex- 
pressions for them in terms of physical quantities for the 
case of a DOPO having a Fabry-Perot cavity of effec- 
tive length L and formed by two identical spherical mir- 
rors with curvature radii R for simplicity. The nonlinear 
crystal is placed at the waist plane of this cavity and 
has refractive index n, second order susceptibility x' 2 ' 
and axial length I (assumed to be much smaller than the 
Rayleigh length of the cavity) . 

For this cavity configuration, the beam radius at the 
waist plane of the resonator is given by [HI 

^ 2= ^V¥ -1 ' j=p,s ' (Ai) 



12 



where Xj is the wavelength of the considered mode inside 
the cavity. 

On the other hand, the parameters which appear in the 
Hamiltonian @ accounting for the external laser pump 
and the nonlinear interaction inside the crystal, namely 
£ p and x, respectively have the following expression: 



X 



2irhc 



(A2) 



h ( c 

Eo \nLXr 



.3/2 



being jj = cTj/2L the cavity decay rate at the considered 
frequency (Tj is the corresponding transmission factor 
through the input mirror), and Pi aS er the power of the 
injected laser. In order to obtain these expressions all 
the transmitted pump power is assumed to be focalized 
inside the transverse dimensions of the nonlinear crystal. 

Through the article, some expressions have been eval- 
uated for concrete system parameters. We have taken as 
typical parameters the following ones 



A p = 400 nm, R 



L = 0.1 m, 



= 1 m, 
1 mm. 



■ 0.1, 
0.01, 



(2) 



2.5, 



Z V ' 



leading to the following model parameters 



w p = 167 /jm, 
7 P = 0.15 ns" 1 



X = 64 s" 1 , 
7 S = 0.015 ns- 1 



having drift vector and diffusion matrix 



and 



V = BB T , 



(B4) 



(B5) 



respectively. 

The parameter v is or 1 depending on whether we in- 
terpret Eqs. (|B1[) as Ito or Stratonovich stochastic equa- 
tions. In the current article any stochastic equation is 
interpreted a la Stratonovich, which allows us using the 
usual rules of calculus. Note that if the extra term in- 
volving derivatives of the noise matrix is zero, Ito and 
Stratonovich interpretations are equivalent, what hap- 
pens for example on Eqs. (fT0|) but not on Eqs. (|22|) . 



APPENDIX C: CONSIDERATIONS ON THE 
LINEARIZATION PROCEDURE 

Following the linearization procedure we explained in 
Section III. A leads not to Eq. ([23]) directly, but to the 
following one 



i {Gh - 2pw ) 9 + b =Ch + gJC (9) C (t) , 



(CI) 



where 



APPENDIX B: CONNECTION BETWEEN 
FOKKER-PLANCK AND LANGEVIN EQUA- 
TIONS 

During the article we make extensive use of the equiva- 
lence between Langevin and Fokker-Planck equations. In 
this appendix we want to briefly review this connection. 

Consider a set of real variables x = (aci, x%, x n ), 
satisfying stochastic Langevin equations 



It 



A(x)+B(x)i,(t), 



(Bl) 



where A (x) and B (x) are a vector and a matrix which 
depend on the variables, and the components of rj (t) are 
real noises satisfying the usual statistical properties (|3T|) . 

The theory of stochastic processes [32j states that the 
stochastic average of any function of the variables / (x) 
can be evaluated as 



(/«> 



stochastic 



ic =y fxp(x)/(x), 



(B2) 



where the probability distribution P (x) satisfies the 
Fokker-Planck equation 



G> r P(r;x) = 



-EM 



(FP) 



2 ^ 

i,3 



P(r;x), 
(B3) 



K{9) =diag (e ie ,e- l9 ,e- l0 ,e l8 ) , 
S = diag (-1,1, 1,-1), 

and the rest of vectors and symbols were defined in the 
corresponding section. Note that the differences between 
this system of equations and the one used in the text (|25|) 
are the matrix JC (8) and the iQh9 term. The latter is of 
order g 2 (as the &'s and 9 are of order g), and hence it 
can be simply removed within the linearized theory. 

Understanding why K. (9) can be removed from the lin- 
earized equations is a little more involved. Projecting 
these equations onto the eigensystem of C ([28]) and defin- 
ing the vector c = col (6, ci, C2, C3) leads to the following 
system of equations (remember that we set cq = 0) 



c = -Ac + gBK (9) rj (r) 



(C2) 



with 



A = 2 diag (0,1,ct - I, a) , 
B = diag(l/2p,j,l,l), 
K(9)=K h3 (9)n 2A (-6), 

being TZij (9) the 2-dimensional rotation matrix of angle 
9 acting on the i—j subspace, and where the components 
of vector r] (t) are real, independent noises satisfying the 
usual statistical properties (13"T]) . Now, we will prove that 
this system, and the same with 1Z (9 = 0) are equivalent 
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within the linearized theory, and hence (|C 1 [) and (|25|) are 
equivalent too. 

To show this, we just write the Fokkcr-Planck equation 
corresponding to this stochastic system (see Appendix 
B), whose drift vector and diffusion matrix are found to 
be 



and 



A = Ac+— col (0,0, 1,0), 



V = g BB 



(C3) 



(C4) 



respectively. Note that in the last equation we have used 
that 1Z (9) is an orthogonal matrix. 

The proof is completed by writing the stochastic sys- 
tem corresponding to this Fokker-Planck equation up to 
the linear order in g, which reads 



-Ac + gBr] (r) 



(C5) 



corresponding to (|C2[) with TZ (6 = 0) as we wanted to 
prove. 

Hence, removing K, (9) and neglecting the iQhO term 
from (|C1[) doesn't change its equivalent Fokker-Planck 
equation within the linearized description, and thus equa- 
tion (|2"5)) must lead to the same predictions as (|C1|) . 



APPENDIX D: LINEAR LANGEVIN EQUA- 
TIONS: SOLUTION, CORRELATION AND SPEC- 
TRUM 



Considering separately the cases r' > r and r' < r, this 
integral is easily carried out, yielding (again the limit 
t ^> A -1 is considered) 



(c(T)c(r')) = -e 



-Al-r'-i 



(D4) 



where we see that this function depends only on the time 
difference \t' — r|, what justifies the name "stationary 
limit" for the t>A _1 approximation. 

Finally, the spectrum of this correlation is found to be 

r+oc _ y2 

C(u)= dfe~ lu}T (c(t)c(t + t)) = 



A 2 + lu 2 ■ 

Hence, particularized to the equations for the projec- 
tions c m ([2"9")h the correlations read 



(d (t) ci (t)) = ' 



5_„-2|r'-x| 

4 



(D6) 



(02 (r)c 2 (r')) = 



9 



4 ((7-1) 



-2(<t-1)|t'-t| 



(03 (r)c 3 (r')>-|^ e - 2 



U T — T 



while the spectra read 



Ci(w) 



9 



4 + w 2 

^2 (W) = — ~ x2 



(D7) 



C 3 (w) 



4(cr- If +0J 2 
9 2 



In this appendix we solve the linear evolution equa- 
tions (|29p of the projections c m , finding their 2-time- 
correlation function and the associated spectrum. These 
equations are of the general type 



c = —Ac + r?? (r) 



(Dl) 



where rj (r) is a real noise satisfying the usual white noise 
statistics pip . A is a real, positive parameter, and T is a 
parameter that might be complex. 

By making the variable change z (t) = c (r) exp (At), 
and considering times larger that A -1 (stationary limit), 
it is straightforward to find the following solution 



c(r) =T / dr 1 r/(T 1 )e A(Tl - r) . 



(D2) 



From the statistical properties of noise P"Tj) we see that 
this solution has zero mean, (c (r)) = 0, and correlation 



(c(t)c(t')) =r 2 e - x ( T+T ') 



(D3) 



r ft 



Jo Jo 



and the stationary limit is reached when t » 0.5, t » 
0.5 (a — l) 2 and r 3> O.Ser -1 respectively. 



APPENDIX E: CALCULATING CORRELA- 
TIONS FOR THE UNDAMPED ORIENTATION 9. 

In this appendix we show how S (n, t 2 ) and C (n, T2) 
defined in (|53[) can be evaluated from the evolution equa- 
tion of 9 (|29a|) . As clearly explained in [33j , this equation 
defines a Wiener or Random Walk process equivalently 
described by the following Fokker-Planck equation: 



d T P(9,r) = ^Dd 2 e P(9,r). 



(El) 



In the same reference, it is also proved that the 2-time 
joint probability associated to this simple diffusion equa- 
tion is 



1 



2D T3-T1 T xi 



2nDy/r 1 (r 2 - n) 



which gives us the probability of passing from orientation 
9\ at time n, to 9 2 at time t 2 (it is assumed that 9 = 
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at the initial time r = 0). Of course, in this expression 
t 2 > T\. From this joint probability the correlations are 
calculated as the integrals 



S(n,T 2 )= d0 1 d0 2 P(6 2 ,T 2 ;6 1 ,T 1 )sme 1 sm0 2 (E3) 



C(n,r 2 ) = [f dOxde^P {0 2 , r 2 ; 6 t , n) cos X cos 6> 2 . 



In order to solve numerically these equations we use the 
semi- implicit algorithm developed in Ref. [3O] - This al- 
gorithm is a finite-differences based method in which the 
total integration time r cn d (the integration is supposed to 
begin always at r = 0) is divided in N segments, creating 
hence a lattice of times {T n } n=0 l N separated by time 
steps At = T cn( \/N. Then, a recursive algorithm starts 
in which the amplitudes at time r„, say /3 n , are found 
from the amplitudes (3 n ~ l at an earlier time r n _i from 



These integrals are easily evaluated, yielding 

S(t u t 2 ) = e-3 D ( Tl+T2 >sinh[£>min(Ti,T 2 )], (E4) 

and 

C{n,T 2 ) = e-5 Z5 ( ri + r2 ) C osh[£>mm(Ti,T 2 )]. (E5) 



APPENDIX F: DETAILS OF THE NUMERICAL 
SIMULATION 

In this last appendix we want to briefly resume the de- 
tails concerning the numerical simulation of the Langevin 
equations (fTS"|) . 

The first important property of these equations is that, 
irrespective of the initial conditions, the amplitudes cor- 
responding to opposite OAM modes become complex- 
conjugate after a short transitory time, i.e., (j3-i, —* 

[/5+i]*V Hence, if the initial conditions are cho- 
sen so that the OAM pairs are complex-conjugate, we can 
be sure that they will remain complex-conjugate during 
the evolution. In particular, we have chosen the above 
threshold stationary solution (j!20|) with 9 — as the ini- 
tial condition. Under these conditions, the 6 Langevin 
equations (|15|) reduce to the following 4 (which we write 
in matrix form): 



with 







A(J3) = 



B(f3) 



, C(r) 



f a-/3 -|/3 +1 | 2 \ 
o-Pt-\Pti\ 2 

5 diag( 0,0, y/fo,\f$ 




(Fl) 



(F2) 



f3 n = (3 r - 



At A 



B[0 



W 



(F3) 



where n is an approximation to the amplitudes at the 
mid-point between r n _i and r n (hence the name "semi- 
implicit" for the algorithm) and the components of W™ 
are independent discrete noises W" with null mean and 
satisfying the correlations 



(Wp, [W?]*> = Ar6 mn 6 jk . 



(F4) 



The mid-point approximation is found from the follow- 
ing iterative algorithm 



AtA (p n ' p - 1 )+B(0 n ' p - 1 



• W n , 
(F5) 

where /3 n, ° =f3 n ~ 1 , being p the iteration index (two iter- 
ations are carried in our simulations), while the discrete 
noises can be simulated at any step as [34| 



(F6) 



with 



r (z, z') = y/— log z cos (2nz') . 



(F7) 



being Zj , z'j , yj and y'j independent random numbers uni- 
formly distributed along the interval [0, 1]. 

This algorithm allows us to simulate one stochastic tra- 
jectory. Then, by repeating it S times, the stochastic av- 
erage of any function can be approximated by the arith- 
metic mean of the values of that function evaluated at 
the different stochastic trajectories. 
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